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We apply the covariant derivative expansion of the Coleman-Weinberg potential to vector-like 
fermion models, matching the UV theory to the relevant dimension -6 operators in the standard 
model effective field theory. The 7 matrix indnced complication in the fermionic covariant derivative 
expansion is stndied in detail, and all the contributing combinations are ennmerated. From this 
analytical result we also provide numerical constraints for a generation of vector-like quarks and 
vector-like leptons. 


I. INTRODUCTION 

The discovery of the 125-GeV Higgs boson with properties consistent with standard model (SM) predictions even¬ 
tually proves the validity of the SM. New physics should reside at a sufficiently high mass scale and/or interact with 
the SM particles through very small couplings, so that it does not affect the SM sector significantly. Thus an effective 
field theory (EFT) of the SM, with all fields appearing in the operators being only the SM ones and all the new 
physics integrated out, would be an ideal tool to study the low energy physics. Today it is well known that going 
beyond the only dimension-5 operator in the seesaw mechanism, up to dimension-6 level such an EFT consists a total 
of 59 independent operators [l[ for one family of fermions as a complete basis, while five more can be added if baryon 
number violation is allowed. 

In the SM EFT collider or other experiment observables can be calculated to the lowest order of the dimension-6 
operators 0, and a global fit @ will determine the Wilson coefficients of the operators in a model independent way. 
On the other hand, if one is interested in a specific model, then translating these model independent constraints 
to model parameters calls for another kind of matching. The most famous example is the Peskin-Takeuchi oblique 
corrections of S and T parameters and the T parameter can even date back to the context of weak charged current 
and neutral current ratio in the form of p parameter. Another example is the Higgs diphoton decay branching ratio 
calculations in the early days of Higgs discovery, hinted by the ATLAS and CMS measured Higgs diphoton channel 
signal strength. In the following we will see how they are equivalent to the Wilson coefficients of the dimension-6 
operators. 

Usually the calculations are done case by case, using the traditional Feynman diagram approach. To the best 
knowledge of the author, the first “systematic” dimension-6 operator calculations are given in [5 Jab Instead of 
the Feynman diagram technique, they use the so called covariant Bin derivative expansion [sl-ll^ (CDF) of 
the Coleman-Weinberg (CW) potential Q, which has the merit of generating all the dimension-6 SM operators 
automatically at one-loop level, including the oblique S and T, the Higgs relevant ones and all the other ones involving 
the Higgs and SM gauge bosons. In this paper we will apply the same technique to a prototype of fermionic modeP, 
namely a mirror vector-like (VL) fermion sector. By definition the mirror VL fermion should contain both doublet 
and two singlets of SU{2)l, but we consider them to be in arbitrary representation (denoted as m, with m = 1,3, 8,6 
etc) of the SU{3)c and the SU{2)l doublet to carry arbitrary U{ 1 )y hypercharge Y. 

Cuided by appendix of the fermionic CDE has extra 7 matrix related contributions. Trace over the 7 matrix 
space immediately picks out combinations of only even number of 7 matrices, leaving odd number ones vanishing. It is 
like another “power counting” in addition to the one of dimension of the operators, and the bosonic CDE corresponds 
only to the case of no 7 matrix. We will see that up to dimension-6 operators there are only two 7 and four 7 matrices 
contribution in extra. 


* e-mail: ran.huo@ipmu.jp 

^ The U parameter defined there corresponds to dimension-8 operator. 
^ For scalar models, further see [T^ . 
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This paper is organized as follows. In Section m we list the diniension-6 operators and the bosonic CDE formula, 
then we generalize it to the fermionic case by working out all possible extra 7 matrix related terms. Then in Section Hill 
we define the VL fermion model and show the main analytical results, as well as some numerical results for the most 
interesting VL lepton and quark model. We discuss the formulism and conclude in Section |TVl At last. Appendix □ 
provides a complete list of contributing combinations, operator by operator in the dimension-6 operator basis and 
term by term in the CDE, in order to facilitate future calculation. 


II. FORMULISM 


Symbol 

Operator expression 

Symbol 

Operator expression 

Symbol 

Operator expression 

Oe 


Ogg 


Ow 


Oh 


Oww 


Ob 


Ot 


Obb 


Ohw 


Or 


Owb 


Ohb 






Od 



TABLE I: Independent CP-even diniension-6 operators composed of only the Higgs and gauge boson fields that are relevant to 
the analysis in this work. Notations of fields and operators are explained in the main text. 


The operator basis we use is listed in Table U which includes only Higgs and electroweak gauge boson fields. For 
operators involving SM fermions generally we don’t need the following technique and a tree level matching [13 will 
suffice. Apparently this is a redundant basis, for example we have 

-{Oww + Owb) — Ow + Ohw = 0 , -^{Obb + Owb) — Ob + Ohb = 0 , (1) 

so that the operator Ohw and Ohb as used in can be switched into Ow, Ob ss well as Oww, Obb and 
Owb- However, the complete list in Table U has an advantage of directly providing a one to one correspondence to 
all dimension-6 operator generated in the CDE calculation. 

Recall that for the bosonic CW potential 


Vcw 


inB 

2 


d‘^P 

{2Tr)^ 


ln(p2 


I""), 


the CDE [6| can be written as a textbook level loop integration 


( 2 ) 
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Here us is the bosonic degree of freedom (DOF) as in the CW potential, 1 for a real bosonic DOF and 2 for a 
complex bosonic DOF. We have assumed dimensional regularization in the loop integration, and in the CDE a Wick 
rotation is performed and subscript “E” indicates Euclidean. V and G generally denote potential and gauge terms 
respectively, and the double prime on V indicates two derivatives with any beyond standard model (BSM) fields. The 
CDE assumes V" = -|- SV" and is some large constant squared mass term which is irrelevant to Higgs vacuum 

expectation value (VEV), SV" on the other hand picks the spacetime dependent part of potential terms. Here not 
only the Higgs VEV but both charged and neutral Higgs components are counted. The SM Higgs sector as well as the 
gauge field can couple between different components of BSM fields, and the SV" should be generalized into matrix, 
with each entry corresponding to coupling term quadratic in BSM component fields. So a trace over the matrix basis 
is needed. ~ generally indicates a Baker-Campbell-Hausdorff (BCH) expansion with covariant derivatives and ^s, for 


® Since the oblique S parameter is given by the dimension-6 Wilson coefficients of S = 4Trv^ (4 cwb + 4- cg), we can see that Ogw 

and Ohb do not contribute to the S parameter. Similarly they will not contribute to the Higgs diphoton branching ratio. 
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(5) 


Here D generally denote a covariant derivative, taking the corresponding gauge field when acting on specific field. 
The generally non-abelian field strength appears as commutator of covariant derivatives of F^^F = ^[Du, Dp\, and 
the Abelian case is easily got. At last /r is a dimension -2 auxiliary number, served for the regularization of the order 
of the expansion. 

In the fermionic case Eq. m becomes 


Vcwo^'-J Hf - I/")- (6) 

On the other hand, the same contribution comes from replacing i) — Vp by —f — Vp. Summing the two logarithms 
and ignoring the —1 inside the logarithm, one gets back to Eq. (12]). Namely the CW potential uniformly describes 
both bosonic and fermionic models. 

We are generalizing the CDE to Eq. ©■ For the eligibility of the DE one also needs a large spacetime independent 
term, so similarly Vp = M + SVp and M is the VL fermion mass. The spacetime dependent 6Vp can only be linear in 

Higgs, up to a dimensionless coupling. In the ln(|i—M—(5y^')-|-ln(— 11 —M—dV^') = \n{{p—M—SVp){;^+M+SVp))-\ -, 

there is always a mixture of dimension -2 and dimension-I terms which corresponds to the bosonic SV” 

sv" = {M,svi;} + isvi;y. (?) 

In this paper {,} and [,] denote anticommutator and commutator respectively. The BCH expansion is obtained 
similarly. 

However, as initiated in the appendix of Q, the presence of 7 matrices introduces further complication, which will 
become our focus from now on. Recall that the CDE always acts on the canonical momentum of Pp + iD^ rather 
than on the kinetic momentum, giving 


e {pp + iDp)F^=J2^[-iD^,[-iD^,[---[-iD-§^,pp + iDp]- 
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The BCH expansion of the first p^ term is acted by the D-^ and is always less in multiplicity of the commutators by 
one than the corresponding action on the second term, so comes the (—-|- l)/(n -|- 2)! factor. The square 
of the expression gives the —G expression defined in Eq. ®, where the first kinetic momentum square term should 
not be counted. On the other hand, the CDE acting on the potential terms is given by the same expression as Eq. Q. 






+ 2^ —r 
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dp dp 


(9) 


Then we only need to use the two equations to calculate 

I + ilp){f + 1$) + iIp){M + 5V^) - (M + + i$) - (M + 5V^ f ) 


-iD-^ 

e ^ 


( 10 ) 


In G we only show terms which can contribute to our dimension-6 operators in Table [U so are the followings. 
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All the relevant terms contributing to the fermionic CDE, including the 7 matrix relevant ones, are generated by 
applying the BCH expansion. 

At first we list the cross term of with M + 6Vp^ which can be regarded as a commutator + 

SVpC'^'^]. Expanding the commutator between terms given by of Eq. © and Eq. (|5]) one by one, 

we get 


±fi = + -{IpDp + D^lp)5V^-§^^ - '-{U)D^D, + Dp$D, + + • • • 




( 11 ) 


The overall sign is up to convention. The first line counts only the commutators involving the first |i in Eq. ®, 
which is proportional to an identity matrix, so that the matrix structure induced commutator vanishes and just the 
^ action on p remains. Note that generally other terms in the expansion has a matrix structure in the BSM particle 
basis, especially the matrix structure of 6Vp and its derivatives do not commute with the gauge generator t“, so the 
second line is nonvanishing. 

Then we consider the two p + ilp terms. We use the identity 7 ^ 7 ^ = ^{ 7 '', 7 '^} + ^[ 7 ^, 7 “"] = g^'^ + ^[ 7 ^) 7 ^], and 
the contraction with the hrst g^'^ always reproduces the bosonic terms. So in addition we have the extra terms of 
1 [ 7 ^, 7 ^] contraction with the two p + iJp expansions of Eq. (|5]), which is 


f2 = -^5[7^7"]l^;.^“ + ••• ■ (12) 

The matrix structure of each term in the expansion is proportional to gauge generator F and identical to two terms in 
the commutator only except the case that the first p is used instead, so the only nonvanishing commutator comes from 
each ^ from the first expansion acting on the Pi, of the second expansion. Eventually very limited terms contribute 
in the CDE. 

In all, the fermionic counterpart of Eq. reads 


^CDE,F = 


np 


du 
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orpE ^ 

(Op)d 2-^ 
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pI 
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(13) 


np = —4 for a Dirac fermion, and the terms to be expanded are in given in Eq. 0 [3 m m [nD. 

Before finishing this section let us describe the way of evaluating the new Ti + r 2 terms. When we identify a 
combination contributing to our target dimension -6 operators, we should first take the trace over the 7 matrix space. 
If one takes the usual 4x4 matrix representation of the 7 matrix as in Q, one have tr( 7 ^ 7 ^) = And this 4x4 

matrix representation induces every term of 5V" + G an internal 4x4 identity matrix in the 7 matrix space, making 
every term multiplied by a factor of 4. On the other hand, we can take the point of view that the overall normalization 
of the CW potential is still determined by the traditional way of counting DOF, such that a Dirac fermion counts 
—4. In this way the internal 7 matrix space is not expanded and we can equivalently take = g^'^. Eq. (fT51) 

is presented in this form. 


III. MODEL AND RESULTS 

Having worked out the general formula, for calculation of a practical model we only need to input the coupling 
matrix 5Vp which is linear in the SM Higgs component fields, and the SM gauge field coupling matrix with the 
matrix structure basically given by the gauge group generator F in corresponding representation. The SM Higgs 
always couples to chiral fermions, and the BSM particles to be integrated out as coupled to each entry of SVp 
should be a Dirac fermion as a requirement of the CW potential, so a chiral projection such as Pp = and 
Pr = ^-^2^ is need. In the calculation performed through Eq. (USD, one should use the symbolic projection algebra 
Pl = Pl,Pr = Pr, PlPr = 0 , Pl -b Pr = 1 repeatedly. 

As a prototype we study one generation of mirror VL fermions, which contains Dirac fermions of one SU (2 )p 
doublet and two singlets. We will neither consider any mixing of the new VL fermions with the SM fermions induced 
by Yukawa interactions, nor bother with any further UV completion issue such as the vacuum stability and the gauge 
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coupling unification. The most general Lagrangian in the interaction basis is 

£ = Ma.ss\ess pure’kinetic terms — Md{udUd + dudi)) — MsuUsus — Msddsds 
- )dsR + yu{u^0L^d^0L){ )usR + yJsLiH~,HS){ 

+h.c.^ 

+ ■^giuo'^^do + doW ud) 

+ \/PT^(uD(i - QuSw)^°^D + rfz>(-5 - QdSw)^°dD + Usi-QuSw)$°Us + ds{-QdSw)^°ds^ 

+ Possible coupling terms with gluon , (14) 


where subscript “D” and “S” indicate the SU{2)l doublet or singlet, and “L” and “R” in Yukawa terms indicate the 
chirality (we changed to Weyl fermion representation for the Yukawa terms). Q« = l ^ Qd = ^ are 

generic electric charge of the upper and lower components of the weak doublet, with Y the hypercharge of the SU (2 )l 
doublet. Note that there are two up types and two down types Yukawa couplings due to the VL nature. 

With the help of chiral projections from the above Lagrangian we can write the fermionic coupling matrix 

0 {yuPR + ruPL)H^ {ydPR + rdPL)H+\ 

Md -{yuPR + yuPL)H~ {ijdPR + ydPL)Ho Qg-, 

-iruPL+yuPR)H+ Msu 0 

{rdPL+ydPR)H^ 0 Msd / 


^ Md 

0 

(ylPL + yuPR)Ho 

V irdPL + ydPR)H- 


The electroweak gauge field coupling matrix, on the other hand, have the upper left 2 by 2 corner filled with SU (2 )l 
gauge field, and the whole diagonal filled with U{1)y gauge field. 

We make the simplification of yu = iu = yu, ijd = Vd = yd and Md = Mgu = Mgd = M. Then the chiral projection 
is trivial, and the VL mass is proportional to identity matrix, always commuting with other matrices. The CDE of 
the CW potential takes the most simplified form. Collecting all the contributions as described in the appendix, we 
get the results 
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Od 
Ogg , 


(16) 


where m = 1,3,8, 6 etc is the representation of the SU{3)c, /(I) = 0, /(3) = —/(8) = —2, /(6) = —| 


® In @ it is pointed out that there are additional universal contributions to the pure gauge dimension-6 operators O 2 S) 

O 2 G OzG defined in the reference. The last two will not affect the electroweak and Higgs physics, while the first three are usually 
small in effect because they are proportional to the SM gauge couplings. For completeness, we also quote the general result here and 
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M [GeV] M [GeV] 


FIG. 1: Constraints on the M-y^ parameter space of the VL lepton (m = 1,Y = —i) model from a global fit to both the 
electroweak oblique corrections and the Higgs data. We have set the neutrino Yukawa = 0. The left plot uses the current 
LEP2 oblique parameters and the current ATLAS and CMS Higgs data. The right plot uses the most aggressive TeraZ result 
for the oblique parameters and the expected FCC-ee Higgs measurement constraints. Green and yellow regions are allowed 
regions at Icr and 2(j level respectively. 


Part of the results can be checked against previous ones. BSM physics contribution to the electroweak oblique 
corrections is related to the diniension-6 operator Wilson coefficient as 

T =- v'^ct , S = 47ru^(4cvi/B + cw + cb) , (18) 

CtEM 

where v = 246 GeV is the SM Higgs VEV. The electroweak oblique corrections are calculated in [H m for a VL 
fermion sector. In the simplification of Vu = Vu = Uu, Vd = = yd and Md = Msu = Msd = M it is not hard to 

expand to order of M~‘^ to get analytical formulas, which agrees with the results given by Eq. (HU) and (|T5)) . It is 
pointed out in [2l|, that the Higgs diphoton decay branching ratio will be enhanced for a VL lepton sector. With 
the help of Table 8, 9, 10, 11 of @ we can connect the above Wilson coefficients to the Higgs precision measurement 
observables, and the charged lepton Yukawa correction to the Higgs diphoton decay amplitude indeed always interferes 
constructively with the SM amplitude, indicating an enhancement of the branching ratio. 


Observable 

fJ^ZZ 

flWW 

y~i~f 

ybb 

fJ'TT 

ATLAS [2^ 

1 44+0-40 
4-^^-0.33 


1.17 ±0.27 

0.52 ± 0.40 

1 40 + 043 

CMS 

1.00 ±0.29 

0.83 ±0.21 

1.12 ±0.24 

0.84 ± 0.44 

0.91 ±0.28 


TABLE H: Signal strengths of various modes, indicated by the subscript in the first column, as measured at the LHC. 

With the dimension-6 operator Wilson coefficients we can fit to the electroweak precision data and the Higgs data, 
to get constraints for the Lagrangian parameter of VL mass M and Yukawa couplings yu and y/'. We do for both 
existing data and future expected measurements. For current measurements we refer to the U = 0 oblique parameter 
measurements of and the ATLAS (CMS) Higgs data with 4.5 (5.1) fb“^ integrated luminosity at ^/s = 7 TeV 


include them in the following fits: 


CD- 


(47r)^ 


(2Y2 + (Y+ l)2 + (y_ l)2)g/2 
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(17) 


For model independent constraints, see [23l. l2‘ 
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FIG. 2: Constraints on the yu-yd parameter space of the VL quark (m = 3, K ~ ^) model from a global fit to both the 
electroweak oblique corrections and the Higgs data. The left plot uses the current LEP2 oblique parameters and the current 
ATLAS and CMS Higgs data. The right plot uses the most aggressive TeraZ result for the oblique parameters and the expected 
FCC-ee Higgs measurement constraints. Green dashed and orange curves are contours at la and 2a levels for different choices 
of M respectively. In the left plot we show contours for M — 200, 400, 600, 800,1000, 1200,1400 GeV and in the right plot we 
show contours for M = 500,1000,1500, 2000, 2500, 3000, 3500, 4000 GeV. 


and 20.3 (19.7) fb ^ integrated luminosity at = 8 TeV , as listed in Table im We do not include tri-gauge 

boson precision measurements in our fitting. For future expected sensitivities, we take the most aggressive oblique 
parameter measurements expected from the TeraZ experiment [2^ and the projected Higgs data from Table 4 of [5^ . 

In Fig. [T] we show the numerical constraints for a VL lepton sector, in which we set the neutrino Yukawa coupling 
yi, = 0. In addition to the oblique S and T parameter constraints, the Higgs diphoton branching ratio calculations 
in [m, now is turned over to a constraint as well. Mostly probing the ye/M direction or so, our results are 
complementary to the direct VL lepton search, which has a mass bound of about 176 GeV [s^. 

In Fig. [2] we show the numerical constraints for a VL quark sector. For small VL mass the constraint curves scale 
more like and indeed the is dominated by Ogg contribution to the gluon fusion in the current LHC data. 

As the VL mass increase the constraint is weaker in the direction of y^ — ydi which is an indication that the oblique 
correction T parameter is becoming more and more dominant in that case. Since every Higgs field is accompanied 
by a Yukawa coupling in the CDE, the Wilson coefficient of a four Higgs operator such as Or will scale quartically 
with the Yukawa coupling, being very sensitive except for the y„ ~ yd direction, in case that the Yukawa couplings 
are large. The direct search constraints for the VL quark masses are 700 ^ 900 GeV [lIHsl. 

IV. CONCLUSION AND DISCUSSION 

The SM EFT is a generalization of the precision measurements observables such as the electroweak oblique correc¬ 
tions and the Higgs precision measurements. At one-loop level the CDE of the GW potential provides a systematic 
tool of matching the UV theory to the SM EFT, all the Wilson coefficients of the dimension-6 operators are calculated 
if we enumerate all the contributing patterns in the CDE. 

The merit of the CDE of GW potential approach includes: 

• Complete: Up to pure Higgs and gauge field operators, all the Wilson coefficients of the dimension-6 operators 
are calculated as long as there is a contribution. On the other hand, it can match UV theory to redundant 
operators, then one can perform field redefinition or use equation of motion to remove redundant operators later 
on, as one wish. 

• Analytical: The CDE do not need to know the mass eigenstates and rotation matrices from interaction eigenstate 
basis to mass eigenstate basis. The complication of mixing is equivalent to the non-commutative nature of the 
matrices of dU" -bG-t-Fi -I-F 2 with each other and with matrices of M~‘^. And the final results of the Wilson 





coefficients in forms of some power of couplings divided by large mass squares, are equivalent to expanding the 
interaction-to-mass-eigenstate mixing matrix to power of 

• Model Independent Formulism: The input is just the BSM particles coupling matrices to the Higgs and SM 
gauge bosons, and all the following calculation are identical for different model. 

Different from the Feynman diagram approach (together with the interaction-to-mass-eigenstate mixing), the CDE 
always gives only the leading terms in expansion of large mass square suppression. The expansion is also precise in 
case that the BSM coupling is small. Even if the current experiments are still exploring a region in which the eligibility 
of the CDE is questionable, as the experiments goes to higher and higher precisions and energies, if no deviation from 
the SM is detected, the eligibility of the CDE will become better and better. 

So far the CDE is only calculated for the case that all the large masses are equal, and the matrix of the large scale 
square is proportional to an identity matrix so that commute with 5V" -|-G-|-ri-|-r 2 . It is not hard to generalize 
it to the case that the large masses are nondegenerate [s^, if we keep track carefully of the order of the generally 
non-commutative matrices. 
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Appendix: List of Contributing Patterns to Each Operator 

Here we collect all the combinations in the CDE of Eq. (fT51) contributing to the dimension -6 operators in Table HI As 
mentioned before the contributing terms can be sorted as 7 matrices irrelevant ones and 7 matrices relevant ones. The 
former group is common in both bosonic theory and fermionic theory, with a correspondence of fermionic {M^SVp} 
and {SVpY' terms to the bosonic counterparts of dimension-1 and dimension-2 parts of 5V". In the following with 
a little bit misuse of terminology we will call the former 7 matrices irrelevant group as “bosonic”, and the latter 7 
matrices relevant terms as “fermionic”. If the theory is indeed bosonic, then one only need to count the pure “bosonic” 
contributions for each operator. 

The Oq is free of derivatives. It is given by all combination of the {5VpY and {M, 5Vp} which make up a total 
dimension of 6 . It can be given by any permutation of three {6Vp)^s, or two {6Vp)'^s together with two {M,6Vp}s, 
or one {SVp)"^ together with four {M,5Vp}s, or six {M,5Vp}s. 

The Oh, Ot and Or form a complete basis for all four Higgs two derivatives terms. In a bosonic theory such terms 
in the CDE come from a permutation of two {SYpYs, or one (SVp)'^ together with two {M, 6Vp}s, or four {M, 6Vp}s. 
A total of two covariant derivatives —should appear at all possible positions; if both act on the same (SVp)'^ 
or {M, SVp}, then an integration by part should be used and one Dp is moved to all the other factors. For fermionic 
theory in addition to the bosonic contribution, there are contributions from permutations of two ilJ)5VpS (first term 
of Eq. (ITT]) 1 together with one {5VpY, and two ilpSVpS together with two {M,6Vp}s. 

Since we keep all the redundant operators and do not use relations such as Eq. m for projection, Ogg, Oww, 
Obb and Owb only receive contributions from the CDE combinations where two Higgs and two gauge field strengths 
directly enter. In a bosonic theory such terms are permutations of a (5VpY together with a 

(fourth term of Eq. ([5])), or two {M,SVp}s together with a jg'^Ft^F^pF^PP -^''. In a fermionic theory additional 
terms are permutations of one {5VpY together with two —jg['jP', 7 '^]F“^,t“s (first term of Eq. (fT^ l. or two {M, 6Vp}s 
together with two — 7 ‘"]F“yt“s, or two ^gjP[SVp, (first term of second line of Eq. (ITTlll. 

The operators Ow, Ob and Ohw, Orb are all two Higgs one gauge field strength and two derivatives operators. 
The integration by parts trick which moves the covariant derivative D to other positions can mix the two groups 
and therefore cannot be used any more. The contributing combinations to Ow and Ob in a bosonic theory are just 
permutations of one ^gpPDDpF^^-^^+ ^gOD^'F^^-^'' (first and second terms of Eq. ([5])) together with one 

—iDp{M,6Vp}-^^ and one {M,6Vp}. In fermionic theory additional contribution comes from permutations of a 
^g^f^lSYp, DpF^pD]-^'^ (third term of second line of Eq. (fTT]) l together with a ilp5Vp, there are also contributions 

from the fermionic two Higgs four derivatives terms. 

The combinations contributing to Orw and Orb are, in bosonic theory a permutation of one igp^FF^^-^ (first 
term in Eq. ®) with two —iDp{M,SVp}-^^. There are a lot of fermionic contributions, from permutations of 
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one together with two ilpSVpS, one ^g^^^[DpSVp, together with one ilpSVp, and 

contributions from the fermionic two Higgs four derivatives terms. 

The bosonic contribution to two Higgs four derivatives operator Od is a permutation of two {M, 5Vp} together with 
four covariant derivatives acting on all possible positions, and the fermionic contribution should be a permutation of 
two ^{IpDp + DplJ))5Vp s (second term of the first line of Eq. (fTTl) '). or one ilpSVp together with one —^{lJ)DpDv + 

DpIpDv +DpDvlJ))5Vp -^^(third term of the first line of Eq. (ITTll b The integration by parts trick can be used to 
always move two derivatives to one Higgs and the other two to the other, in accord with the form of Od- Note that the 
four covariant derivatives are eventually contracted by two metric g^'^s, and the desired Od corresponds to a specific 
order of the two D^s and two DyS. In order to get to the desired order one need to commute Dp and Dy. In bosonic 
theory all the gauge field strength from the commutator of [Dp, Dy] = —igFI^^^D has been collected and resummed 
by the expansion of G, so the commutator should not be counted. But in fermionic theory the resummation of gauge 
field strength terms has not been done in the 7 matrix related term, not all of them are collected by the G or r 2 - 
The [Dp, Dy] = igFpyD in switching orders among the four derivatives should be counted, which makes additional 
contribution to the Ow, Ob and Orw, Orb- In particular, organizing terms to the form of two derivatives acting 
on one Higgs and two on the other, there can be three different orders: directly {DpDf^H^){DyD''H) = Od, as well 
as [DpDyF[^){D^D'^F[) and {DpDyF[^){D''D^Fl). {DpDyF[^){D^^D''FI) will have a contribution to Ow, Ob as well 
as Orw, Orb, and [DpDyH'^){D''D^^H) will have a contribution only to Orw, Orb- To evaluate the contribution 
a transformation based on Jacobi identity is found useful. 
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